Charge dynamics in half-filled Hubbard chains with finite on-site interaction 
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We study the charge dynamic structure factor of the one-dimensional Hubbard model with finite 
on-site repulsion U at half filling. Numerical results from the time-dependent density matrix renor- 
malization group are analyzed by comparison with the exact spectrum of the model. The evolution 
of the line shape as a function of U is explained in terms of a relative transfer of spectral weight 
between the two-holon continuum that dominates in the limit U oo and a subset of the two-holon- 
two-spinon continuum that reconstructs the electron-hole continuum in the limit U ^ 0. Power-law 
singularities along boundary lines of the spectrum are described by effective impurity models that 
are explicitly invariant under spin and jy-spin SU{2) rotations. The Mott-Hubbard metal-insulator 
transition is reflected in a discontinuous change of the exponents of edge singularities at ?7 = 0. The 
sharp feature observed in the spectrum for momenta near the zone boundary is attributed to a Van 
Hove singularity that persists as a consequence of integrability. 
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I. INTRODUCTION 

Since its proposal^ the Hubbard model has become 
a paradigm in the field of strongly correlated electron 
systems. It is the simplest model that accounts for the 
metal-insulator transition on a half-filled lattice when the 
on-site electron-electron repulsion U is strong enough. R 
is still debated whether the model in two spatial dimen- 
sions or some variation of it contains the mechanism for 
high temperature superconductivity at finite doping. 

Theoretically, much more is known about the model on 
a one-dimensional (ID) lattice.^ In this case, it is possi- 
ble to calculate the exact spectrum and eigenfunctions by 
Bethe ansatz (BA) 4 Two remarkable properties revealed 
by the exact solution are the existence of fractional ex- 
citations that carry separate spin and charge quantum 
numbers and the opening of the Mott-Hubbard gap at 
half filling for arbitrarily small U > Q. 

Recently, there has been renewed interest in dynami- 
cal properties of ID models. One motivation for this is 
that questions about features of the excitation spectrum 
of ID systems, such as the persistence of spin-charge 
separation at high energies, have become relevant with 
the improvement in the resolution of momentum-resolved 
experiments In addition, ultracold atoms trapped in 
optical lattices have emerged as a new means to study co- 
herent dynamics of ID models, including integrable ones 
which are not realizable in condensed matter systems.^ 

At the same time, significant progress has been 
achieved in developing analytical^""^^ and numerical^^"— 
techniques to study dynamical correlation functions in 
the high energy regime where conventional Luttinger liq- 
uid theory^®'^^ does not apply. Analytically, it is possi- 
ble to compute exponents of power-law singularities that 
develop near thresholds of the spectrum of dynamical 



correlation functions at arbitrarily high energies. For the 
metallic phase of the Hubbard model, i.e. away from half 
filling, the calculation of finite-energy dynamical correla- 
tion functions was pioneered by the pseudofermion dy- 
namical theory. This theory is based on the BA solu- 
tion of the model and has been applied to calculate, for 
instance, the optical conductivity and the one-electron 
spectral function of ID conductors;^"— In another ap- 
proach, exponents of high energy singularities can be in- 
vestigated using effective field theories that treat high en- 
ergy modes as impurities, defined in momentum space, 
which can scatter off low energy excitations (see Ref. [22 
for a review). This approach, combined with the BA so- 
lution, has also been applied to calculate edge exponents 
for the spectral function of the Hubbard model away from 
half filhng 

In this work we are interested in finite energy dynam- 
ical correlation functions for the Hubbard model at half- 
filling. Clearly, the edge exponents of the Mott insulating 
phase should differ from those of the metallic phase stud- 
ied in Refs. [ToIITtI due to the finite charge gap. In fact, 
the result should be simpler due to the higher symmetry 
of the model at half filling. At half filling the Hubbard 
model has a hidden r/-spin SU (2) symmetry that rotates 
between doubly occupied sites and empty sites.— In the 
same sense that spin SU{2) invariance fixes the expo- 
nents of spin correlation functions at high energies,^* it 
should be possible to use the continuous symmetry in 
the charge sector to constrain the exponents in charge 
dynamics. 

Particularly, we shall focus on the charge dynamic 
structure factor (DSF) S{q,uj) at zero temperature. The 
DSF is known analytically only in two limits. In the low 
energy limit, which requires that the Mott-Hubbard gap 
be small, dynamical correlation functions can be calcu- 



2 



lated using form factors for the integrable sine-Gordon 
modelii^ Within this field theory approach, a universal 
square-root cusp is found at the edge of the relativistic 
spectrum of massive charge solitons. In the strong cou- 
pling limit, one can take advantage of the factorization of 
the wave function into a noninteracting charge sector and 
a spin sector described by the ID Heisenberg model;^ 
A square-root cusp is again found at the lower thresh- 
old of the two-holon continuum, stemming from matrix 
element for noninteracting spinless fermionsi^ In both 
low energy and strong coupling limits, no features are 
predicted at the branch line of the spin excitations here 
called spinous. This is in contrast to the behavior of the 
one-electron spectral function, which has sharp features 
near both charge and spin branch lines for any value of 

The DSF has also been studied numerically; ^^'^^ most 
recently using the dynamical density matrix renormaliza- 
tion group (DDMRG).^'* Most of the numerical work has 
focused on the regime of large U, which is appropriate to 
describe strong Mott insulators such as Sr2Cu03»^ 

We note that the DSF has strikingly different line 
shapes in the weak and strong coupling limits. For non- 
interacting electrons, U — 0, the DSF can be calculated 
exactly and corresponds to the density of states of an 
electron- hole pair. For U oo, the spectral weight is 
assigned to a two-holon continuum, with negligible con- 
tribution from spinous. ™ 

The purpose of this paper is to investigate the charge 
DSF for the Hubbard model at half filling for arbitrary 
values of g, at finite U . We construct a picture for the 
intermediate U regime by combining information about 
the exact spectrum from BA, an effective field theory for 
edge singularities at high energies and numerical results 
from the time-dependent density matrix renormalization 
group (tDMRG). We start in Sec. [IT] by discussing the 
exact support of the DSF in terms of elementary charge 
and spin excitations using known results from the BA so- 
lution. Our main results can be found in sections IIIII and 
IIVI In Sec. IIIII we present the effective field theory that 
incorporates the spin and 77-spin SU{2) symmetries ex- 
plicitly and allow us to determine the exponents of power- 
law singularities at the edges of the exact spectrum of the 
DSF. In Sec. |IV]we present the tDMRG results for cer- 
tain values of U and analyze them by comparison with 
the field theory combined with the exact spectrum from 
BA. In addition, we discuss the U dependence of the line 
shape, interpolating between weak and strong coupling 
limits. Finally, Sec. IVl contains the conclusions. 

Our results are relevant for the charge DSF of fermionic 
atoms in a ID optical lattice with on-site atomic re- 
pulsion described by the integrable Hubbard model. In 
the context of cold atoms, the charge DSF is probed by 
Bragg spectroscopyi^Si^ The results are also useful as 
an approximation to condensed matter systems where 
the intcgrability-breaking perturbations to the Hubbard 
model, such as the nearest neighbour interaction in the 
extended Hubbard model, are small. In this context 



the DSF has served to interpret electron energy loss 
spectroscopy^ and inelastic x-ray experiments. ^^'^ Since 
the Shiba transformationiS maps the charge DSF for 
[/ > to the spin DSF ioi U < 0, our results also apply to 
the spin DSF of the spin-gapped phase for the attractive 
Hubbard model. 



II. 



MODEL, SYMMETRIES AND EXACT 
SUPPORT OF THE CHARGE DSF 



Model 



We consider the ID Hubbard model 



H 



L 

E 



h.c. 



U 



(1) 

Here Cj = (cj,^, Cj,i) is a two-component spinor represent- 



ing electrons with spin a - 
and L is the system size, 
site j is denoted as Uj = 



:t, i at site J, Uj^rj 
The number of electrons at 



We focus on half 



filling (rij) = 1. We have set the hopping amplitude to 
1, not to be confused with the real time variable t. The 
ID Hubbard model is integrable. The exact spectrum for 
all values of f7, density and magnetization is provided by 
the BA solution.^ 

At half filling and zero magnetic field, the Hubbard 
model has an explicit spin SU{2) symmetry and a less 
obvious charge Ty-spin SU{2) symmetry.^ The generators 
of spin rotations are the components of the usual spin 
operator S = J^j '^j('''/2)Cj, where t is the vector of Pauli 
matrices. The generators of 77-spin rotations are 



I "^{nj.t + nj.i - 1) = ^ 77|, 



E(-i) 



(2) 
(3) 
(4) 



such that the local operators obey the algebra [rijTrjp] = 

±5jjirj^ and [77^,77^^] = 26jj>rij. Notice that 77^ is pro- 
portional to the fiuctuation of the local charge density 
operator: rij = 1 + 2r]j. The transverse components rj^^y 
of the 7/-spin vector ry = {r]'" , rj^ , rj^ ) can be defined by 
= 77^ ± irjy . 

While the spin and 7;-spin symmetries account for an 
5(9(4) = [SU{2) X SU{2)]/Z2 symmetry, the global sym- 
metry of the Hubbard model was recently found to be 
larger and given by [SO{4) x f/(l)]/^2 = [SU{2) x 
SU{2) X C/(1)]/Z|.4i In addition, the ID model has an 
infinite number of local conserved quantities associated 
with integrability. 
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B. Charge structure factor at half fiUing 

The charge DSF is defined as the Fourier transforni of 
the density-density correlation function 

S{q,io) - ^ E \{GS\n,W)\H{io - + Egs) 

v^GS 
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die-*(7?|(t)7?o^(0)), (5) 



\GS) is the ground state and 



where n, 

\v) is an excited state with energy E^. Since S{q,uj) — 
S{—q,uj), in the following we set q > without loss of 
generality. The DSF obeys the sum rule 

r ^ ^(9,^) = 2(n,,^n,,;) ^ 22), (6) 

with the density of doubly occupied sites given exactly 



1 + cosh(a;[//2) " 



(7) 



For U — 0, we have T) — ^ 1/4; for U — > oo, the integrated 
spectral weight vanishes as 2) — >■ In2(2t/L/)^. 

The ?7-spin symmetry can be used to relate the DSF 
at half filling to the correlation function for the pairing 
operators, which create doubly occupied or empty sites. 
The ground state for an even number of sites is unique 
and is a singlet of both spin and ?7-spin rotations (quan- 
tum numbers S = = rj = i]^ = 0). Eq. ([5]) can be 
rewritten as 

S{q,uj)^87TL J2 \{GS\rfo\iy)\''5g^pJ{uj - E, + Eos), 

u^GS 

(8) 

where is the lattice momentum of the eigenstate \v). 
By employing, for instance, the unitary transformation 
that rotates the 77-spin vector by 7r/2 about the y axis, 
U = e~'^^^ , we can rewrite the matrix element 

{GSWo\v) = {GS\U^Ur^^„U^U\v) = {GS\r^^„\v') , (9) 

where \v') = Ulv) is also an eigenstate of H with energy 
E^i = E^, but with momentum P^i = P^, + tt. The 
momentum shift follows from the fact that the lattice 
translation operator anticommutes with 77^.— We then 
have 

= SnL J2 \{GS\vSW}\Hg,p^+^ 

u^GS 

xd{uj - E, + Egs) (10) 



4 g-'"('?+^)3 



the charge DSF tensor 

/oo 
dte-*(f (t)r)g(O)), (11) 
J -co 

where a,b = x,y^z and f]^ = r]j, fjj'^ = (— 1)^?]^'^. In 
this notation, S{q,uj) = S^^{q,uj). The 77-spin SU{2) 
symmetry implies 

S{q,iu) = 5+-(g,c.)/2 

/OO 

Therefore, up to the shift of total momentum by tt, the 
line shape of the charge DSF is identical to that of the 
correlation function for pairing operators rf^ . We can 
also write 

5(g,w) = -Ee-^« / dte-*(i7,(t). 170(0)). (13) 



For later reference, we mention that for [/ = the 
charge DSF in Eq. ([5]) reduces to the density of states 
for excitations with a single electron-hole pair 



'S'o(g,w) 



Ae{uj - uj^{q))e{uj+{q) - Lo) 



(14) 



and likewise for the correlation function for Thus 
S{q^uj) can be viewed as the longitudinal component of 



where uj^{q) = 2 sing' and uj+{q) = 4sin((7/2) are the 
lower and upper thresholds of the electron-hole contin- 
uum, respectively, and 9{uj) is the Heaviside step func- 
tion. Up to a factor of 2, this is the same result as for 
spinless fermions at half filling. The free electron DSF 
has a step discontinuity at the lower edge and a square- 
root divergence at the upper edge, which stems from the 
Van Hove singularity of an electron and a hole with the 
same velocity. 



Elementary excitations in the Bethe ansatz 
solution 



In this subsection we review some BA results for the 
exact spectrum which will be useful for comparison with 
numerical results in section IV. 

According to the BA solutionji^ the eigenstates of the 
ID Hubbard model can be constructed from elementary 
charge, 77-spin and spin excitations. In the half filling 
case, it suffices to consider two branches of excitations, 
one in the charge sector, which we call holons, and one 
in the spin sector, which we call spinous. (For the rela- 
tion between the holons and spinous used here and the 
notation used e.g. in Refs. I4l1l43l . see Appendix A.) 

In the thermodynamic limit holons and spinous have 
dispersion relations edpc) and es{ps), respectively, where 
the dressed momenta pc.s and dressed energies ec,s are 
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given by [see Ref. S and Appendix A; here we follow 
the notation in Eq. (7.8) of Ref. 2] 



Pc{k) 
edk) 



2 " 



TT 

2 



f°° Jo(w) sin(w sin fc) 

- 2 / aUJ ; TTTT-- — . 

°°^^ yo(a;)sin(a;A) 
Q w cosh(wC//4) ' 



2cosfc + U/2 

Ji(cj) cos(tj sin k) 



-4 I duj 

lo ^(1 



Ji{lo) cos(a;A) 



/o Cij cosh(a;t//4) 



(15) 
(16) 

(17) 
(18) 



Here k and A are the charge quasimomentum and spin 
rapidity, respectively. For any value of U > 0, the spin 
dispersion is gapless at the spinon Fermi points ps — O^n 
and the charge dispersion has minimum energy at pc — 
— 7r/2, with a gap given by 



duj 



1 



sinh(27rtj/t/)' 



(19) 



Analytic expressions for the holon and spinon dispersions 
can be obtained in the limits [/ — > oo: 



£c(j>c) 



U 

2" 

4 

U 



2 sinpc 



41n2 
U 



(1 + cos2pc 



smps 



and in the limit [/ — !■ 0: 

£c{Pc) ~ 
£s{Ps) ~ 



2pc 



2smps, 



(20) 
(21) 

(22) 
(23) 



with an exponentially small charge gap A w 
(4/7r)V[7e~^'^/^. We shall also be interested in the ve- 
locity of holons and spinous, defined as 



UciPc) 



dec 
dpc' 



Us{ps 



dps 



(24) 



The dispersion relations of holons and spinous given by 
Eqs. (dtl-IHD are illustrated in Fig. [T] 



D. Boundary lines in the exact spectrum of S{q, uj) 

Even though the exact spectrum and wave functions 
of the ID Hubbard model are known, it has not been 
possible to calculate the DSF directly from the BA solu- 
tion. The difficulty is in computing the matrix elements 
in Eq. ([5]) for significantly large chains. Unfortunately, 
unlike the Heisenberg model, there are so far no determi- 
nant formulas'*^ or vertex operator approach^ to com- 
pute form factors for the Hubbard model. 

Nonetheless, we can use the BA equations to com- 
pute the exact support of the DSF. It follows from the 




Figure 1: (Color online.) Exact dispersion relations of ele- 
mentary charge and spin excitations for U = 4.9. Here the 
excitations are represented by particles in the empty holon 
(c) and spinon (s) bands. The holon band is gapped with 
minimum energy at pc — —n/2. The spinon band is gapless 

at Ps = 0, TT. 



Wigner-Eckart theorem that the excited states that con- 
tribute to S{q,LL!) in Eq. ([S]) must carry quantum num- 
bers S — = (spin singlets) and ?/ = 1, f]^ = (7/-spin 
triplets). This selects states with 2m holons, m > 1, and 
2n spinous, n > 0. Since the excited states must contain 
at least two holons and the holon dispersion is gapped, 
the DSF vanishes for uj < 2 A. 

The simplest excited states, in the sense of lowest num- 
ber of elementary excitations, that contribute to S{q, w) 
are two-holon states (m = 1, n = 0). For rj = 1, i]^ = 0, 
the excitations with m — 1^ n — have total momentum 
P and energy E given by^ 

P ^Pc.l+Pc,2+T^, E = £c{Pc,l) + £c{Pc,2), (25) 

where and Pc.2 are the dressed momenta of the indi- 
vidual holons as in Eq. ([T5|) . The next simplest excited 
states that contribute to S{q,io) contain two spinous in 
addition to the two holons (m ~ n — 1). For 77 = 1, 
77^ = excitations with m ~ n ~ 1, we have 

P = Pea +Pc,2 +Ps:2 + 7r, (26) 
E = edPci) + SciPca) + £s{Ps,i) + £siPs,2), (27) 

where Ps^i and Ps,2 are the momenta of the two spinous 
as in Eq. p^ . We expect these two classes of states 
to give the leading contributions to the spectral weight 
of S{q,uj) for all values of U, based on the observa- 
tion that, analogously, the leading contribution to the 
half-filling one-electron excitations stem from one-holon- 
one-spinon excited states.''^ Indeed, figure 2 of Ref. \^ 
presents the contributions of different states to the one- 
electron-addition sum rule for half filling. Interestingly, 
the higher-order contributions are most important at 
U ~ 4, yet they account only for about 0.005 of the 
one-electron-addition spectral weight. Consistently, it is 
expected that the higher-order contributions associated 
here mainly with to = 1, n = 2 and m — 2,n — states 
are again very small and maximum at U ~ A. 
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Figure 2: (Color online.) Support of the charge DSF for U = 
4.9. Energies are measured from the Mott-Hubbard gap 2A. 
Some special lines of the spectrum are shown. The two-holon 
continuum is bounded by (-U2c{<l) and 0^2^(5), but the spectral 
weight is nonzero everywhere above the lower line cj^^2s il) ^^^"^ 
extends to arbitrarily high energies. For < g < go, where 
go is the momentum of the point indicated by a diamond, the 
lower edge of the two-holon continuum (jJ2cil) defined by 
two holons with the same momentum (— tt + g)/2. For go < 
q < n, the energy of two holons with momentum (— tt + g)/2 
follows the dashed line, but this is no longer the lower edge 
of the two-holon spectrum. Instead, a;^;,(g) is defined by two 
holons with different momenta but equal velocities. For < 
g < g<i, where g< is the momentum of the point indicated 
by the left-pointing triangle, the DSF vanishes below a;^j,(g). 
For g< < g < g[> = TT — g<], the absolute lower threshold is 
ijJ2c2sil)' defined by an excitation with two holons with the 
same momentum, one spinon at the Fermi surface and another 
spinon below the Fermi surface that has the same velocity as 
the holons. For g[> < g < tt, the line (jJ2c2s{i) defined by two 
holons with the same momentum and two spinous at opposite 
Fermi points. 



Fig. [2] illustrates the exact support of the DSF. It also 
indicates special boundary lines in the m = 1, n = and 
m = n ^ 1 continua which will be important to construct 
the effective field theory for edge singularities in Sec. IIIII 
as well as to analyze the tDMRG data in Sec. IIVI 

We now discuss the most important boundary lines in 
the spectrum of S{q^uj) based on simple kinematics. For 
large U, we expect the spectral weight of S{q,u!) to be 
confined inside the two-holon continuumi-'^ The upper 
threshold of the two-holon continuum uJ2c(.l) given by 
two holons with the same momentum (tt -|- q)/2, < q < 
TT. In the strong coupling theory for [/ — >■ cxi,^'^ in which 
limit the holons have a free-fermion cosine dispersion, 
the lower threshold of the two-holon continuum is given 
by two holons with the same momentum pc,i = Pc,2 = 
{—TT+q)/2 for all < q < TT. However, for any finite U the 
holon dispersion deviates from the cosine function such 
that the curvature of the dispersion (absolute value of 
inverse effective mass) is smaller near the minimum of the 
band than near the maximum. As a result, for values of q 
near the zone boundary the two-holon excitation with the 



lowest energy has holons with different momenta pc,i = 
q—T^—Pc,2 7^ Pc,2 (mod 27r), but such that they propagate 
with the same velocity, Uc{pc,i) — Uc{pc.2)- Starting from 
Pc,i = Pc,2 = — 7r/2 and increasing the holon momenta, 
the values ofpc,i andpc,2 that define uj2^{q) split off at the 
inflection point of the exact holon dispersion. Thus there 
is a value of q^^iU), given by twice the momentum of the 
inflection point (plus or minus tt as in Eq. (I25p ). where 
the nature of the lower threshold changes. In the limit 
[/ ^ 00, Eq. ^ yields q^ ~ tt - 161n2/[/ + 0(C/-2). 
Using the exact holon dispersion in Eqs. (jTS]) and (|17p we 
find that decreases monotonically with U and ^ 
in the limit U ^ 0. 



The lower edge of the two holon continuum is not the 
absolute lower threshold of the support of S{q, oj) for gen- 
eral q. Starting from q = and moving along the line 
uj2ci<l)j ^ value of q is reached at which the velocity of the 
holons with momentum (— 7r + g)/2 becomes equal to the 
spin velocity at the spinon Fermi surface. The value of 
q = q<{U) where this happens is given by the condition 
Mc((— TT + <z<j)/2) = Vs = Us{0) and is represented by a 
left-pointing triangle in Fig. [21 For q > q^, it is possible 
to lower the energy by transferring momentum to a pair 
of spinous. For q^ < q < n — q^ = q^^ the lower edge of 
the two- holon- two-spinon continuum, denoted i^2c2sil)-: 
has two holons with pc,i = Pc.2 < {^tt + <7)/2, one 
spinon at the Fermi point with p^ i = and another 
spinon with momentum ps,2 = q ~ 2pc.i such that the 
velocity of the latter equals the velocity of the holons, 
UsiPs.2) — Uc{Pc,i)- For q^ < q < TT^ the lower edge has 
the two spinous pinned at opposite Fermi points while 
the holons carry the same momentum q/2. 



The line '^2c2s{l) actually the absolute lower edge 
of the support for q^^ < q < i:. Adding more holons to 
the excited state can only increase the energy due to the 
charge gap. Furthermore, we find numerically that the 
spinon band has no inflection points away from the Fermi 
surface. In this case the minimum energy for 2n spinous 
at fixed total momentum is obtained by placing 2n — 1 
spinous at the Fermi surface and one spinon carrying the 
remaining momentum, giving the same minimum energy 
as for two spinous only. Notice that the i^2c2s{^) li'^^ 
not the same as the spinon mass shell, in contrast with 
the lower edge for the metallic phaseJii^S 



Finally, we note that in the limit U ^ Q the line 
^2c2s il) becomes the lower edge of the electron- hole con- 
tinuum, (^2c2sil) ~^ 2 sing, whereas the lower edge of 
the two-holon continuum becomes the upper edge of the 
electron-hole continuum, uj2cil) ~^ 4sin(g/2). As U — >■ 0, 
we expect that all the spectral weight of S{q, uj) becomes 
confined between UJ2c2sil) ^'^"^ '-^2cil) order to recover 
the free electron result. 
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III. SU{2) INVARIANT IMPURITY MODEL 
FOR EDGE SINGULARITIES 

In this section we work out the field theory methods 
that allow us to describe power-law singularities of dy- 
namical correlation functions at high energies. The gen- 
eral method relies on effective impurity models to treat 
the high energy modes. This approach has been applied 
to other models and is explained in detail in Ref. [l^. 
Here our goal is to extend these methods to incorporate 
the spin and 77-spin SU{2) symmetries of the Hubbard 
model at half filling explicitly in the effective impurity 
models. The main idea is to define vector currents for 
the high energy modes, in analogy with the low energy 
SU{2) currents used in the Sugawara representation of 
the spin part of the Luttinger model4^ 

A. Low energy theory 

Before dealing with high energy singularities, we re- 
view standard results obtained by bosonization of the 
Hubbard model in the low energy limit, The starting 
point is to linearize the electron dispersion for U = 
about the right (R) and left (L) Fermi points for the 
two spin channels a =t,i- In the continuum limit, the 
fermionic field is expanded in the form 

c,,a ^ ^ e'^"/Vi?,^(x) + c-'^^/^^lA^)- (28) 

Bosonization maps the fermionic fields to 

^cA^) - i^c.,.e-'^^-"("\ (29) 

for a — L,R = , where Fa^^ are Klein factors. 
The chiral bosonic fields satisfy [(pa,u{x),dx'ipa' ^v'ix')] = 

1016 a a' ^v^v' 

5{x-x'). Charge and spin bosons are defined 
as the linear combinations 

'Pa,c{^) ^ VPaA^) + VaA^)]/^^ (30) 
ipaAA ^ VPaA^) ~ VaAA]/^- (31) 

The long wavelength part of the spin and r^-spin density 
operators can be expressed in terms of the chiral spin and 
charge bosons as 

S, ^ S{x)^3rAx)+3lAx), (32) 
J7j ■q{x) ^ 3rA^) + ^lA^)^ (33) 

where Jq.i/ with v = c^s are SU{2) charge and spin cur- 
rents with components 

JAix) = a9,^„,,(x)/V4^, (34) 
J±Jx) = e^^Vi^'P.,Ax)/2n. (35) 

These SU{2) currents obey the k — 1 Kac-Moody 
algebra4^ We remark that the long wavelength parts of 
S{x) and ri{x) do not mix charge and spin bosons, but 
the staggered parts omitted in Eqs. (15^ and ([55)) do.-^ 



In the low-energy limit, spin-charge separation holds 
in the strong sense that spin and charge excitations 
are decoupled. The bosonized version of the Hubbard 
model in Eq. ^ yields the Hamiltonian density 7i{x) — 

E.=c,J^^°H2^)+'5H.(x)] with 

^(0) ^ ^(j|^^+j2^^)^ (36) 

SV-u = -2T:VuK3Ft.,v ■ (37) 

The terms u'^A are quadratic in the bosonic fields and 
can be recognized as the Luttinger model for charge 
and spin collective modes written in manifestly SU{2) x 
SU (2) invariant form. The parameters Vc and Vs are the 
charge and spin velocities, respectively. For U ^ 1, we 
have Vc ^ 2 + U/2Tr and Vs ^ 2 — U/2Tr. The terms SHi, 
are perturbations that mix R and L currents and are not 
quadratic in the bosonic fields. For J7 ^ 1, Ac ~ — C//27r 
and As ~ U/27T. Although the bare coupling constants 
A^ are small for U ^ 1, these perturbations flow under 
the renormalization group with /3 function 

^^-A + 0{Xl), (38) 

where d£ — \dA\/A with A the high energy cutoff. For 
U > OAs is marginally irrelevant and the spin spectrum 
is gapless. On the other hand, Ac is marginally relevant 
and gives rise to a charge gap. The gap A ~ e-^/\^<:\ 
is exponentially small at small U, in agreement with the 
BA solution (c.f. below Eq. The charge sector 

can then be described using the sine-Gordon modelj^ 
whose elementary excitations are solitons with a massive 
relativistic dispersion e{q) — A {vcQY + Note the 
roles of spin and charge bosons are exchanged if we invert 
the sign of [/, as follows from the Shiba transformation!^ 
The critical theory of the spin sector is the k — 1 SU (2) 
Wess-Zumino-Witten (WZW) model.'*^ In the more ele- 
gant notation of non-Abelian bosonization, operators can 
be written in terms of the 2x2 unitary matrix field g{x) 
of the WZW model, 

AxA = ^aA^A® 9r{xA. (39) 

where x± = Vgt ± x and the tensor product notation 
means gij = gj^ i9Rj with i,j = 1,2. The chiral spinor 
fields gL and gR have conformal dimensions (-1,0) and 
(0, j)A respectively, and can be represented in Abelian 
bosonization notation as 

5"(^)= [ e»v^^„,.W j ' a = L,R. (40) 

Under a spin rotation represented by a unitary 2x2 ma- 
trix U, the chiral spinors transform as ga^. g'ai — 
Uijga.j- Due to conformal invariance, the spin SU{2) 
symmetry is enlarged to a chiral SU{2)l x SU{2)r sym- 
metry. In terms of the matrix field, the spin currents are 
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given hy^ 

3L,s = ^TTid+gg^T), J^,, - -^Tr(gta_gr), (41) 

where d± — d/dx±. The theory for the low energy sec- 
tor of the Hubbard model is equivalent to that of the 
Heisenberg spin chain, the only distinction being in the 
spin velocity Vs, which depends on U. 



B. Edge singularities at high energies: imposing 
spin SU{2) invariance in spin correlation functions 




Figure 3: (Color online.) (a) "Deep hole" particle-hole excita- 
tion that gives the lower edge of the longitudinal spin DSF. 
In the effective field theory, spinons are interacting spinless 
fermions with hole states for — vr < p < and particle states 
for < p < TT. (b) Particle-hole excitation with high energy 
particle. Due to spin inversion symmetry, the dispersion is 
particle-hole symmetric and (a) and (b) are degenerate. 



Although low energy theories based on the linear dis- 
persion approximation yield reliable results for thermo- 
dynamic quantities, in general they fail to predict the cor- 
rect edge singularities of dynamic correlation functions^ 
For this purpose it is important to take into account for- 
mally irrelevant perturbations that break the Lorentz in- 
variance of the fixed point Hamiltonian. Nonlinear Lut- 
tinger liquid theory makes progress by refermionizing 
the elementary excitations. For spin-1/2 models, this 
means defining spinless fermions associated with holon 
and spinon bands that have a finite curvature about the 
Fermi points. ^^'^^ 

The idea behind the effective impurity models for edge 
singularities is the same for all dynamic correlation func- 
tions. Essentially, it involves defining high energy sub- 
bands within the dispersion of elementary excitations, in 
addition to the chiral low energy modesi^i The single- 
particle states used to define the high energy sub-bands 
depend on the momentum and energy of interest for the 
dynamic response function. In order to motivate the ap- 
plication of the SU (2) invariant effective field theory for 
edge singularities, let us turn for the moment to the case 
of spin correlation functions, for which more is known 
concerning the implications of SU{2) invariance J^ii^ We 
will show that the proposed definition of a high energy 
impurity spinor in Eq. (|48p below recovers known results. 



1. Lower edge of the two-spinon continuum 

For the half-filled Hubbard model with J7 > 0, the 
spectrum of spin correlation functions is gapless. The ef- 
fective theory for edge singularities of the spin DSF has 
been worked out for the XXZ model, which only has 
U{\) symmetry for general anisotropy parameter but in- 
cludes the SU{2) symmetric Heisenberg point. In the 
spinless fermion language, the spin excitations are de- 
scribed by particles and holes in an interacting band (see 
Fig. [3]). The longitudinal spin DSF is defined as 

/oo 
dte"^'{S^{t)S^m. (42) 
-OO 



We can also consider the transverse spin DSF 

/oo 
dte^-\S+{t)So{Q)). (43) 



Spin SU{2) invariance at zero magnetic field implies 
SI%q,0j) = S+-iq,io)/2. 

The lower edge of the support of S^^{q, w) corresponds 
to the lower threshold of the two-spinon continuum and 
is described as a "deep hole" excitation with a hole with 
momentum p = —q below the Fermi point and a particle 
exactly at the Fermi point. The energy of this excitation 
is equal to the spinon mass shell Ss{q) > 0. Since at zero 
magnetic field the spin band is particle-hole symmetric)^ 
the excitation with a hole at the Fermi point and a par- 
ticle ai, p = q above the Fermi point is degenerate with 
the deep hole excitation. 

The edge singularity in this case is described by a q- 
dependent effective model which, besides the low energy 
states near the Fermi points, contains impurity sub-bands 
associated with the deep hole or the high energy parti- 
cle. The spin DSF shows a singularity above the spinon 
mass shell, S^^{q,LL>) ^ Slu'^, with Suj = uj — es{q)- The 
lower edge exponent /i is determined by the scaling di- 
mension of the operator that creates the particle-hole ex- 
citations after performing a unitary transformation that 
decouples the impurity modes from the bosonized Fermi 
surface modes. For details, see Ref. [lH. After this uni- 
tary transformation, up to irrelevant operators, the ef- 
fective Hamiltonian density assumes the noninteracting 
form n = + nf\ where 

y.^^ = dlies - iusdx)ds + d\{ea - iusda:)ds. (45) 

Here, da{x) and da[x) are field operators that annihilate 
a high energy spinon particle and a deep spinon hole, 
respectively, and < Vg is the velocity of both impurity 
sub-bands. The high energy sub-bands are defined with 
momenta centred at ±g and have momentum cutoff A, 
with Us A <C 1 (see Fig. [31). The ground state is a vacuum 
of dg and dg ■ After the unitary transformation, the spin 
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operator that is applied to the ground state is of the form 

-(I^(x)e*^''«'^«'=("^)e*^''^'^^'=(="\ (46) 

where the relative minus sign between the two terms 
comes from ordering the Klein factors of the sub-bands 
(recall creates a hole). For the U{1) symmetric 
model, the parameters jr_l can be related to exact phase 
shifts.-'^ In the case of SU{2) symmetry, these parame- 
ters can be fixed by the condition that longitudinal and 
transverse spin correlations have the same exponents.™ 
This condition implies = 1/^/2 and 7l = and the z 
component of the spin operator reduces to 

S'{x) cx 4(a;)e-'^'^«'=(^) - dl{x)e'-^'^"^^'^'=\ (47) 

The dimension- 1/4 vertex operators in Eq. (|47l) can be 
recognized as the components of the chiral spinor gji in 
Eq. (|40l) . This observation motivates regarding d and d 
as the components of a high energy spin impurity spinor 



dj,{x) 
ds{x) 



(48) 



which must transform under spin SU(2) rotations as 
Ds^i —?' D'^ ^ — UijDs.j- With this definition, the particle- 
hole degree of freedom of the impurity is interpreted as 
an effective pseudospin 1/2. The operator in Eq. (|T7)) 
can be rewritten in the compact form 



S'{x)^Dl{x)T'gR{x). 



(49) 



In fact, the equivalence of longitudinal and transverse 
correlation functions follows from the correlation func- 
tions of the spin SU (2) vector operator 



S{x) cx Dl{x)TgR{x). 



(50) 



The transverse components in Eq. (j50p also agree with 
known results J^ii^i^ 

The free Hamiltonian in Eqs. (jHJ and (j45|) can be 
rewritten in the SU{2) invariant form 



n 



27rw, 



(51) 



In this effective model for the lower edge singularity, the 
states in the Hilbert space are constrained to have ei- 
ther zero (ground state) or one impurity (excited states), 
Nd,s = J dxDl(x)D^[x) = 0,1. There is no essential 
distinction between the two high energy sub-bands since 
the transverse components of the total spin vector 



S = J dx [Jr,s{x) + Jl,s( 



{x) + Dl{x){T/2)D,ix)] (52) 



generates rotations of deep holes into high energy parti- 
cles. The time ordered propagator for the free Dg field 
reads 

{TD^/x,t)Dl^{0,0)) = S,,,9{t)e-''^'S{x~u,t), (53) 





Figure 4: (Color online.) (a) Particle-hole excitation that 
gives the upper edge of the two-spinon continuum in the lon- 
gitudinal spin DSF. (b) 5^ = -1-1 excitation corresponding 
to the creation of two spinous above the Fermi points. Due 
to spin inversion symmetry, the latter is degenerate with the 
S" — excitation in (a). 



where Dg^i = ds and Ds,2 ~ dg. The correlation func- 
tions for the chiral spinors are given by the standard con- 
formal field theory result 

{9lA^,t)9LA0,0)} cx 5.,,(x+)-i/2 (54) 
{9RA^.t)9RA^M « S,,,{x^)-'/\ (55) 

Using these expressions, we can calculate the edge ex- 
ponent from Sl^{q,u}) J dx J dt e'"^* {B{x,t)B^O,0)) 
with B^^x) oc Dl{x)T'gR{x). This gives n = -1/2, the 
same as the result for the Heisenberg modeli^ 

In order to connect with the methods developed for 
U{1) symmetric models, Hamiltonian (|44| must be inter- 
preted as the effective model after the unitary transfor- 
mation that decouples the mobile impurity. However, a 
different approach could be to write down Eq. (jSip di- 
rectly based only on SU (2) symmetry. In this case, in 
addition to the terms in Eq. (|5T|) we would be led to 
write down the marginal operator 



6T-Lrld = -2TrVs{KR3R^s + hlJl,s) ■ dItD^ 



(56) 



where kr^l are dimensionless coupling constants. The 
longitudinal part of this operator amounts to a density- 
density interaction between the impurity and the Fermi 
surface modes. The full operator SHrld is equivalent to 
a two-channel Kondo coupling, which appears naturally 
in the problem of a mobile spin- 1/2 impurity coupled to 
a ID electron gas.— In Appendix |B] we show that the 
Kfl,L operators are marginally irrelevant for k_r,l > 
(equivalent to ferromagnetic Kondo coupling). Although 
we are not able to derive the bare coupling constants 
starting from the Hubbard model for general [/, we shall 
assume that hr^l are positive for [/ > because oth- 
erwise we would not recover the known results for the 
Heisenberg model. Moreover, it is known that the finite 
size spectrum for excited states of the Hubbard model 
that contain high energy holes in the spin band fits the 
"shifted" conformal field theory form^ii suggesting that 
the marginal operator should be irrelevant for any fi- 
nite U . With the asymptotic decoupling of the impurity 
spinon, the symmetry of the effective model (|51l) becomes 
SU{2)l X SU{2)r X SU{2)d. 
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2. Upper edge of the two-spinon continuum 

The SU{2) invariant effective theory can aiso be ap- 
plied to tlie upper edge of the two-spinon continuum, 
where it is known that the spin DSF for the Heisen- 
berg model has another power law singularityji^ In this 
case the threshold is given by a particle with momentum 
q/2 and a hole with momentum — q/2, as shown in Fig. 
3^. In this case the excited state has two impurities. 
The particle and hole states form the components of a 
single impurity spinor Dg as given by Eq. ()48p . Thus 
Nd,s — J dxD\{x)Dg{x) — 2 for excited states. We in- 
troduce the time reversal conjugated spinor 



-dlix) 
dUx) 



(57) 



which transforms like Dg under spin rotations. The ex- 
cited state that describes the upper threshold of the two- 
spinon continuum is created by acting on the ground 
state with the operator B^'^{x) oc DI{x)t^ Dl{x), where 
the high energy particle and high energy hole in the final 
state must be treated as distinguishable particles, as in a 
two-body problem. SU{2) symmetry dictates that the 
effective impurity model is of the form 



2'KV, 



(58) 



Here we have included the parabolic term in the disper- 
sion of the impurities, with effective mass < 0. 

The marginal part of the operator in Eq. (|58p acts 
on the excited state as a density-density interaction be- 
tween the two impurities. For J7 > 0, we expect Va < 
as obtained for the Heisenberg model,— implying an at- 
tractive interaction between particle and hole. The Vs 
interaction turns out to be crucial for the upper edge 
singularity Sa^{q,uj) ^ 6uj'^, with Sut — 2es{q/2) —ui. For 
Vs ^ 0, the density of states diverges as (5w — > due 
to the Van Hove singularity for particle and hole with 
equal velocities. However, for any Vs ^ 0, the solution 
of the two-body problem shows that the matrix elements 
are strongly affected by resonant scattering and turn the 
divergence into a square-root cusp with /i = +1/2. The 
effect is analogous to a ID exciton problem for particles 
with negative mass, hence no particle-hole bound state 
above the continuum for Vs < 0. 

But what we have described is the interpretation of 
the singularity in the longitudinal spin DSF. An alter- 
native route to determine the edge exponent would be 
to rely on the spin SU (2) symmetry and consider the 
transverse spin DSF. In this case, instead of a particle- 
hole pair, the excited state has either two particles with 
momentum q/2 (for S~~^ {q,uj)) or two holes with mo- 
mentum —q/2 (for Sg {q,oj)) (see Fig. Uja). The ex- 
cited state with — -1-1 is created by the operator 



case of the transverse components we need to intro- 
duce the point splitting because the operator creates two 
spinless fermions with approximately the same momen- 
tum. Thus the leading term has higher scaling dimen- 
sion than the longitudinal component B'^^ . On the other 
hand, for spinless fermions the Vs interaction is irrele- 
vant — the s-wave scattering amplitude vanishes — and 
can be neglected in the effective Hamiltonian. Remark- 
ably, we encounter the same exponent fi — 1/2 due to 
matrix elements for free spinless fermions with vanishing 
relative momentum. '^'^ This can be verified by calculating 
the propagator for the pairing field dldxdl}^ Therefore, 
SU{2) symmetry tells us that the upper edge exponent 
can be interpreted as due to either strong interactions in 
the excitonic pair or statistics of free spinless fermions. 



C. Edge singularities at high energies: imposing 
r;-spin 5(7(2) invariance in the charge DSF at half 
filling 

We now turn to edge singularities in S{q, w), which in- 
volve the creation of high energy holons. Within the field 
theory approach, we represent the charge excitations as 
holes in a completely filled band or particles in an empty 
band, with Mott-Hubbard gap 2A. We will borrow the 
nomenclature often adopted in the literature and refer to 
these bands as the lower Hubbard band and the upper 
Hubbard band, respectively. Since there are no Fermi 
points in this case, the holon band only contributes with 
impurity sub-bands to the effective model. By analogy 
with Ds in Eq. (|48p , we define the charge impurity spinor 
for given high energy holon sub-bands as 



(59) 



such that creates a particle in the upper Hubbard band 
and dl creates a hole in the lower Hubbard band. The 
ground state is a vacuum of Dc- Due to 7y-spin SU{2) 
symmetry, explicit in Eq. p3p . the effective Hamilto- 
nians as well as the operators that create high-energy 
excitations in the field theory must be written in terms 
of the charge impurity spinor. The generator of ry-spin 
rotations is represented by 



dxDl{x){T/2)D^{x). 



(60) 



{x) cx DI{x + |)T+£)*(a 



d\dxd^^- In the 



We are now in a position to compute the exponents for 
the thresholds of the charge DSF in Fig. [5] 



1. Boundary line co2cil) /°'' 9 < 9« 

Consider first the lower edge of the two-holon con- 
tinuum for momentum in the range g < go, such that 
^2cil) — 2£c(— 7r/2 -|- q/2). The effective model in 
this case has two charge impurities in the excited state, 
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g/2 q/2 




E 




Figure 5: (Color online.) (a) Particle-hole excitation that 
gives the lower edge of the two-holon continuum in the charge 
DSF for q < q^. In the effective field theory, holons are spin- 
less fermions with a gap between the lower Hubbard band and 
the upper Hubbard band, (b) rf = excitation that adds 
two particles to the upper Hubbard band. Due to particle-hole 
symmetry, the latter is degenerate with the rf =Q excitation 
in (a). 



iV^^^ = J dxDl{x)D^{x) = 2. The 77^ = state corre- 
sponds to a hole in the lower Hubbard band and a parti- 
cle in the upper Hubbard band, as illustrated in Fig. [S^. 
The particle and hole are the components of the same Dc 
spinor and the situation is analogous to the upper edge 
of the two-spinon continuum in the spin DSF. Due to the 
ry-spin SU{2) symmetry in Eq. (|13p. the edge exponent 
can also be calculated from the excited state of two parti- 
cles created in the same sub-band (Fig. [Sb). The vector 
operator that creates these ?7-spin triplet excitations is 
B^x) (X Dl{x + ^)tD*{x - §). The effective Hamilto- 
nian density consistent with SU{2) x SU{2) symmetry 
reads 



n = DI i Ec- iucdx 



92 
2m. 



D^ + V,iDlDj^ (61) 



Due to symmetry, there is no coupling between holons 
and low-energy spinous at the level of marginal operators. 
Since rric > 0, we expect Vc > for absence of a particle- 
hole bound state below the threshold. It follows that 
the edge singularity is of the form S{q, u) ^ 5lo^ with 



5uj 



„{q) and /i — 1/2. A similar conclusion can 



be reached for the singularity at the upper edge of the 
two-holon continuum ^^((j) for all values of q. We note 
that ?7-spin rotations mix states with 77^ = 0, ±1, but 
the total momentum of the r]^ = state differs from the 
momentum of the ry^ = ±1 states by tt. This is consistent 
with the spectrum from the BAi^ 



2. Boundary line cj2c2s il) q< < q < q> 

For q^ < q < q^. — IT ~ q^, the lower edge of the sup- 
port of S{q, uj) has one low energy spinon and one impu- 
rity spinon in addition to the two holons. The operator 
that creates this two-holon-two-spinon excitation must 
be constructed using one low energy chiral spinor, one Ds 
spinor and two Dc spinors. Furthermore, selection rules 
impose that the operator is a vector of ?7-spin rotation 
and a scalar of spin rotation. These conditions naturally 
lead to B^x) oc DI{x+^)tD*{x- ^)Dl{x)gR{x) as the 



operator with the lowest scaling dimension. Besides the 
sum of Eqs. ([?T|) and (|61l) with Uc — Us — u, the effective 
Hamiltonian contains the symmetry allowed interaction 
between the spin impurity and the charge impurities 



(62) 



The parameter Vcs could in principle be related to the ex- 
act phase shift in the nontrivial S matrix between a high 
energy holon and a high energy spinon. We then need to 
compute the propagator for three impurities that move 
with the same velocity, interact among themselves but 
are decoupled from the low energy modes. It is easiest to 
discuss the rj^ — +1 excitation instead of the 77^ = one, 
trading the interactions between distinguishable charge 
hole and charge particle by the problem of noninteract- 
ing holons which are indistinguishable fermions. Simple 
power counting in the correlation function for (x) (the 
calculation is detailed in appendix [C]) yields the edge 
singularity S{q,Lu) ^ Suj^ with 5uj = u — ^■2c2s(l) ^^^"^ 
H = 3/2. 



3. Boundary line i^2c2s il) qo < q < iv 

For q^ < q < TT, the lower edge of the support has two 
spinous at opposite Fermi points. Thus we are looking for 
a spin scalar operator that involves the low energy modes 
only. The momentum tt scalar operator of the WZW 
model is the trace of the matrix field Tr[g(a;)], which 
has scaling dimension 1/2. The operator that creates the 
excitation in this case is then B^ix) oc Dl{x+ ^)tD^{x — 
|)Tr[(ji(x)]. Again, we find the edge exponent /i = 3/2. 



4- Boundary line u!2^{q) for q^ < q < iv 

Finally, let us discuss the lower edge of the two-holon 
continuum for q^^ < q < tt. In this case the 77^ = excited 
state has a hole in lower Hubbard band and a particle in 
the upper Hubbard band that move with the same veloc- 
ity, but are not associated with the same charge impurity 
spinor. We denote the spinor for the holon with momen- 
tum below the inflection point of the holon dispersion (see 
Fig. [1]) by Dc and the spinor for the holon above the in- 
flection point by Dc- The occupation of the impurity sub- 
bands in the excited state is Nd.c ~ J dx D\{x)D^{x) — 1 
and Nd.c — J dx Dl{x)D^{x) = 1. The vector operator 
in this case reads B^x) oc DI{x)tD*{x). The effective 
Hamiltonian density with marginal operators allowed by 
symmetry is 



n 



Dl 



iucdx 



2m, 



iUcdx 



2rhr 



(63) 
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where and are the Coulomb and exchange in- 
teractions between the distinguishable impurities. (For 
models ((58|) and (|61|) with a single impurity spinor, these 
two interactions are equivalent.) 

The model in Eq. (|63p is again similar to a ID exciton 
problem. There is a Van Hove singularity in the density 
of states when the relative momentum between Dc and 
Dc holons approaches zero. We expect this divergence to 
be removed for arbitrarily weak final-state interactions. 
There is a priori no reason why V^' and should be 
zero or even small at finite U . Depending on the sign of 
the effective scattering amplitude, a bound state can be 
formed below the continuum, which is in fact observed 
numerically for the extended Hubbard modeL^^ 

However, in appendix |D] we show that the existence 
of nontrivial conservation laws in the Hubbard model 
requires — Vf' — exactly. Remarkably, the inte- 
grability of the model implies that impurity holons asso- 
ciated with different Dc spinors do not scatter off each 
other. 

We stress that the vanishing of and does not 
follow from 77-spin SU{2) symmetry alone. This is rea- 
sonable because it is possible to generate infinitely many 
models with the same symmetry that are not integrable, 
for instance by adding finite range ry-spin exchange in- 
teractions J2j j' Jj-j'Vj ' Vj' to the Hubbard model. For 
non-integrable models, we generically expect the forma- 
tion of two-holon bound states^^ below cj2cil > Qo) — as 
well as the broadening of any power-law singularity that 
is not protected by kinematics. 

When we set = ~ 0, the propagator of B^^(a;) cx 
DI{x)tDI{x) factorizes into free propagators for Dc and 
Dc impurities. The Van Hove singularity of the density 
of states persists in the DSF as S{q, lu) ^ Su'-'- with 6uj = 
(jj — Ec — £c and /i = —1/2. This is the only divergent 
edge singularity in the charge DSF and only appears at 
finite energies and finite U. 

The results for the boundary lines discussed in this 
section are summarized in Table HI The exponents for 
the lines UJ2c{q < 9o) and CLi^^(q) are consistent with the 
large-C/ results of Ref. 33. The exponent for the line 
^2c2s(l> < <Z < tt) agrees with the low energy result 
obtained assuming Vc = Vg in Ref. Isil . Our results 



Table I: Predictions of the SU{2) invariant effective impurity 
models for the charge DSF of the Hubbard model. The bound- 
ary lines considered here are illustrated in Fig. [5] As the 
frequency approaches the boundary lines, Su — >■ 0, the DSF 
behaves like S{q,u}) oc / da; / e*"* (B(a;, t) ■ (0, 0)) ~ Slj''. 



Boundary hne 


Vector operator 


Edge exponent /i 




DlrDt 


1/2 




DItD* 


-1/2 


^2ci<l) 


DItDI 


1/2 


^2c2si<l-i <q<q>) 


DItDIDUr 


3/2 


'^2c2s{l > 1>) 


DlTDlTT[g] 


3/2 



show that these exponents hold at finite U and away 
from the low energy limit. The exponents for the lines 
a; 2c (go < (7 < tt) and (^2c2s('Z4 < 9 < 9t>) could not be ob- 
tained by either large-[/ or low energy approximations. 
Notice that the exponents predicted by the SU{2) in- 
variant impurity models are all half- integers, in contrast 
with the continuously varying exponents of the metallic 
phasciiiSiiZ 



IV. NUMERICAL RESULTS 
A. Methods 

We have used the tDMRG method to compute the 
real time density-density correlation function G'(j, t) — 
(nj(t)no(O)) for Hubbard chains with open boundary 
conditions and lengths up to 200 sites. The method starts 
with a traditional DMRG calculation, ^2*53, obtaining the 
ground state \GS) of the finite chain. The single site 
operator r^g for a central site is applied to the ground 
state, and then this state is evolved in real time, ob- 
taining \ip{t)). The original ground state is retained in 
matrix product state (MPS) form, so that the tDMRG 
need only target |V'(0)- each time step we measure 
('7j (^)'7o(0)) by measuring the off-diagonal MPS overlaps 
{GS\rij\ip{t)) for all sites j. A single run provides results 
for all frequency and momenta by Fourier transforming 
over time and space (i.e. j). 

The time evolution operator is written as a product of 
exact nearest-neighbor bond exponentials, as in a famil- 
iar Suzuki- Trotter breakup. Recently Kirino, Fujii, and 
Ueda have reported excellent performance with a partic- 
ular fourth order breakup, in which every bond operator 
is applied in every half-sweep, but in reverse order for 
every other half-sweep.^- We have also found that this 
method gives very small finite-time-step error and ap- 
pears to be superior to other breakups for high accuracy 
calculations. 

The main limitation of the tDMRG method is on the 
maximum time reached by the simulation, due to the 
growth of entanglement with running time. Typically, we 
have reached ^max ~ 20 in units of inverse hopping, keep- 
ing a maximum of m = 2500 — 4000 states. We find that 
the entanglement grows more rapidly for smaller values 
of U and this prevents us from studying U < 1. The spa- 
tial Fourier transform is done first, and no windowing is 
required since within the maximum time reached, the sig- 
nal which is propagating within \ip{t)) has not yet hit the 
edges of the system. Thus the resolution in momentum is 
not limited by the system size. Windowing is necessary in 
the time Fourier transform, but the frequency resolution 
would be poor if we fit the window within i,nax- Instead, 
we extrapolate the time signal using linear prediction, al- 
lowing the use of a larger window The resulting line 
shapes for the charge DSF do not have any analytic in- 
put. A conservative estimate for the frequency resolution 
of these line shapes is given by l/^max- This resolution 
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Figure 6: (Color online.) Charge DSF of the Hubbard model 
at half filling calculated by tDMRG as function of momentum 
q and energy a; for [/ = 4.9. 



could be substantially improved by using analytic results 
for the edge singularities of the DSF to help extrapolate 
the DMRG data to much longer times. 



B. tDMRG results for ^(g, cj) 

We now analyze tDMRG results for {7 = 1, C/ = 2, 
and U = 4.9, obtained without any analytic input, by 
comparing with the predictions of the field theory in Sec. 
mil combined with the exact spectrum from the BA. 

First we discuss the result for S{q,uj) for U — 4.9 
shown in Fig. [51 The exact support of the DSF in this 
case is illustrated in Fig. [2l notice, however, that the 
energies in Fig. [2] are shifted by the Mott-Hubbard 2A 
while the energies in Fig. [6] are not. The spectral weight 
distribution in Fig. [5] is consistent with the strong cou- 
pling picture'^^ in the sense that the spectral weight is 
rather small below the lower threshold of the two-holon 
continuum. However, for values of q near the zoire bound- 
ary it is already visible that the onset of the spectral 
weight occurs below the lower edge of the two-holon coir- 
tinuum. As discussed in Sec. Ull the main contributioir 
to this weight is due to excitations with two spinous in 
addition to two holons and the support of S{q, w) extends 
down to the line i-^2c2sil)- 

Another featured observed in the tDMRG results for 
U = 4.9 is a sharp asymmetric peak above the lower edge 
of the two-holon continuum for q near the zone boundary. 
This effect is predicted by the theory in Sec. IIIII as a 
change in the exponent of the edge singularity from /i ~ 
1/2 for uj2c{q < Qo) to ^ — —1/2 for u}2ci<l > 9o)- Using 
the exact holon dispersion for U — 4.9, we obtain go « 
0.447r. Fig. [7] shows constant-g cuts of S{q,uj) for q = 
O.Stt < go and q — O.Gtt > go- The arrows indicate the 
threshold energies predicted by the BA. 




Figure 7: (Color online.) Line shapes of S{q, uj) calculated by 
tDMRG for U = 4.9 and two values of q. The arrows indicate 
the exact edges of the spectrum predicted by BA. The field 
theory predicts a square-root cusp at w = oj^ for q — O.Stt 
but a square-root divergence at w = for q = 0.67r. 



In order to confirm the existence of two regimes for the 
ijJ2c{q) edge exponent, we have analyzed the time decay of 
the momentum depeirdent correlation function G{q,t) — 
e~"'-'G(j, t). We assume an asymptotic power-law 
decay of G(g, t) and fit the real part G(g, t) in the time 
range 7 < t < 20 to the formula 



Re G(g, t) = Aq cosiWqt + (f>q)r 



(64) 



with Aq,Wq,(f)q,riq sls frcc parameters. Since S'(g,aj) is 
given by a time- frequency Fourier transform of G{q,t)^ 
the exponent rjq in G(g, t) is related to the exponent fi in 
5(g, uj) ^ 5lj^ by ?7q = 1 + /i for the smallest /i among the 
boundary lines. The fitting to Eq. (|64|) should work best 
in the range go < g < tt, in which we predict a square- 
root divergence in S'(g, w) which strongly dominates the 
long time behavior of G(g, t). 

The time decay of Re G(g, t) is illustrated in Fig. [S] 
The energies and exponents obtained by fitting the nu- 
merical results to Eq. (|64p are shown in Fig. [51 We 
first note that the frequencies extracted from the tDMRG 
data are in excellent agreement with the exact result from 
the BA. In fact, the tDMRG are slightly shifted to higher 
energies as expected from the error due to the finite Trot- 
ter step. 

Furthermore, the results for the exponent in Fig. [5| 
clearly show 7]q w 1/2 for g near the zone boundary. 
This supports the existence of a square-root divergence 
in S'(g, uj) which corresponds to the Van Hove singular- 
ity predicted by the theory in section [HI] as due to the 
absence of scattering between distinguishable impurities 
in the integrable model. We note that the existence 
of a bound-state below the continuum would lead to a 
non-decaying contribution to G(g,t), which is not ob- 
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Figure 8: (Color online.) Real part of the density-density 
correlation function G{q,t) — X]j s"''"* ('^i(*)'^o(0)) for U — 
4.9. Symbols represent tDMRG data and solid lines are fits 
to the power-law decay form in Eq. (|64[l for 7 < t < 20. The 
data for q = n/4 has been rescaled by a factor of 10. 
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Figure 10: (Color online.) Charge DSF for U = 2. 
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Figure 9: (Color online.) Frequencies (circles) and exponents 
(diamonds) obtained by fitting the tDMRG results for G{q, t) 
for U = 4.9 to Eq. (|64|l . The dashed line represents the exact 
lower edge of the two-holon continuum. 



served. On the other hand, the error in the numerical 
value of the exponent increases with decreasing as the 
energy window of validity of the square-root divergence in 
S{q,uj) decreases, which implies that longer times would 
be needed in order to observe the asymptotic behavior of 
G{q,t). Nonetheless, Fig. [5] suggests that the exponent 
is significantly larger below q = ~ QAAn. Recall that 
the prediction of the effective impurity model is = 1/2 
for q < q^, which gives rjq = 3/2. 

Let us now discuss the result for U — 2 shown in Fig. 
[TUl For this smaller value of U, we see that a larger frac- 
tion of the spectral weight is located below the two-holon 
continuum. The lower edge of the support agrees with 
the exact line U!2c2sil) ^'^^ U = 2. We can quantify the 
distribution of spectral weight by computing an average 



frequency Lj{q) from the first- moment sum rule as 
dLJLuS{q,Lu) _ idtG{q,t = 0) 



u{q) 



J^dujS{q,u;) Giq,t = 0) 



(65) 



The expression on the right hand side of Eq. (|65|) is 
directly provided by the tDMRG from the short time be- 
havior of G{q,t). For U = 4.9 the average frequency 
Lj(q) is always above a;^^((7). In contrast, for U = 2 we 
find that a7(g) < w^^(g) for q > 0.647r. The difference 



„(q) — uj{q) increases as q 



Therefore, it appears 



that the small U behavior, characterized by all the spec- 
tral weight lying below the two-holon continuum, is ap- 
proached more rapidly for larger values of q. 

The transfer of spectral weight to below the two-holon 
continuum as U decreases is confirmed by the result for 



U — I showir in Fig. [TTJ In this case the lines 



and uj2ciq) are already very close to the lower and up- 
per thresholds of the electron- hole continuum for U = 0, 
respectively. However, there is still significant spectral 
weight in the two-holon continuum. 

The results in Figs. [TUl and [TT] reveal that S{q, uj) has a 
rounded peak below the lower edge of the two-holon con- 
tinuum. The peak is more clearly seen in Fig. 1121 which 
shows the line shape for g = tt for J7 = 1 and U = 2. Par- 
ticularly in the case U = 1 the peak is very narrow and 
the spectral weight is rapidly suppressed below the oir- 
set of the two-holon contribution. We can also see that, 
although a large fraction of the total spectral weight is 
associated with two-holon-two-spinon states, the singu- 
larity above 1^2^(9) persists. 



C. General picture for the DSF at finite U 

111 light of the analytic results in Sec. IIIH the line 
shapes in Fig. [T2]suggest a scenario for the U depeirdence 
of S{q, uj). When combined with the exact spectrum from 
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Figure 11: (Color online.) Charge DSF for U — 1. 





Figure 13: (Color online.) Spectral weight for states that 
contribute to S{q,Lj) for a 10-site Hubbard ring for q = tt [(a) 
and (b)] and g = 47r/5 [(c) and (d)]. The center of each line 
represents the energy of an excited state as a function of 
1/U or U. The intensity is proportional to the matrix element 
I {GS\nq\i')\. In (a) and (c), the open triangles on the left-hand 
side mark the energies corresponding to free spinless fermions 
in the large U limit. In (b) and (d), the open triangles mark 
the energies of one-particle-hole excitations in the U = case, 
and the filled triangles the energies of the two-particle-hole 
excitations. 



Figure 12: (Color online.) Line shape of S{q,uj) calculated 
by tDMRG for q = n and U — 1,2. The exact free electron 
result for (7 = is also shown for comparison. 



the BA, the SU{2) invariant effective field theory does 
not predict any divergence below the lower edge of the 
two-holon continuum. However, the free electron result 
in Eq. (|14p exhibits a Van Hove singularity from below 
the upper threshold of the electron-hole continuum. We 
interpret Fig. [12] as indication that the free electron line 
shape is recovered as the peak below a;^^((7), which is 
rounded for any finite U, becomes narrower as J7 — 0. 
Only at [/ = does the Van Hove singularity develop 
at what is then the upper threshold of the electron-hole 
continuum. 

Moreover, for any finite U and fixed q > Qo (recall that 
go — >■ for J7 — 0) the square-root divergence above 
w^^(g) is always present. However, the spectral weight in 
the two-holon continuum vanishes for J7 — > 0. The total 
spectral weight of S{q, lu) is not conserved as U varies (see 
Eq. ([6])), but in relative terms the weight is transferred 
from the two-holon continuum for [/ — > oo to the subset 
of the two-holon-two-spinon continuum that lies below 



the lower edge of the two-holon continuum for [/ — >■ 0. 

The subset of the two-holon-two-spinon continuum 
that dominates S{q, w) and reconstructs the electron-hole 
continuum in the limit ?7 — >■ can be obtained from the 
heuristic rule that the holons are constrained to the min- 
imum of the holon band (momentum pc — — 7r/2 in Fig. 
[1]), where the Mott- Hubbard gap closes for [7 = 0, while 
the spinons are free to move along the spinon band. We 
conjecture that for J7 — the matrix elements for the 
charge density operator in Eq. ([5]), which are not known 
except for small chains, select excited states with two 
holons and two spinons according to this rule. 



D. Lanczos results for small systems 

In order to provide further evidence for the above sce- 
nario, we have calculated S{q, uj) for a 10-site half-filled 
chain with periodic boundary conditions by exact diag- 
onalization based on the Lanczos method. Figs. [T5k 
and 113b illustrates the energies and matrix element for 
all eigenstates of the Hamiltonian with total momentum 
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q = -IT. The important point is that for this small system 
there is only one state that gives a large contribution to 
S{q — TT, uj) in both limits of large U and small U. This 
is the state that has energy equal to 4 at C/ = 0, which 
corresponds to the maximum energy for an electron-hole 
excitation with q = n. 

By solving the Lieb-Wu equations^ for system size 
L = 10, we have computed the exact energies of two- 
holon states and identified that the state that evolves into 
the upper edge of the electron-hole continuum at [/ = 
is the lowest energy two-holon excitation. All states 
with energy lower than the latter involve excitations in 
the spinon band. This observation is consistent with the 
proposed scenario for the U dependence of S{q, uj) since 
it shows that the state that defines the lower edge of the 
two-holon continuum and carries a large spectral weight 
splits off from the continuum below it for arbitrarily small 
U . In the thermodynamic limit we expect that this be- 
havior corresponds to the disappearing of the Van Hove 
singularity below the upper edge of the electron-hole con- 
tinuum and the formation of another Van Hove singular- 
ity above the lower edge of the two-holon continuum once 
we turn on the interaction. 

We have also calculated the matrix elements for ex- 
citations with momentum q = 47r/5 for the chain with 
L = 10 (Figs. and[T3Jl). Interestingly, for < g < tt 
there is a level crossing as a function of U where the 
spectral weight associated with the lowest energy two- 
holon state changes abruptly. This is a manifestation 
in the small system of the change in the nature of the 
lower edge of the two-holon continuum from = 1/2 
to ^ — —1/2. The value of U where the level crossing 
happens is given by the condition qo{U) — q aX fixed g, 
where qo(U) is twice the value of the momentum at the 
inflection point of the single holon dispersion. Indeed, 
in Fig. [T3t we see that the weight in the lowest energy 
two-holon state is larger on the small U side of the level 
crossing {1/U > 1/8), which corresponds to the regime 
where we expect a square-root divergence above ^2^(5) 
in the thermodynamic limit. 



V. CONCLUSION 

In summary, we have studied the charge dynamic 
structure factor S{q, oj) of the Mott insulating phase of 
the ID Hubbard model at finite U, based on a combina- 
tion of Bethe ansatz, field theory and tDMRG. 

We used the BA solution to discuss the exact spec- 
trum of excitations that contribute to S{q,uj), without 
low energy or strong coupling approximations. Unlike the 
metallic phase, the lower edge of the support of S{q, u) 
is not given by the spinon mass shell, but by either the 
lower edge of the two-holon continuum or the lower edge 
of the two-holon-two-spinon continuum that has three 
particles (two holons and one spinon) at finite energies 
with the same velocity. In addition, an important differ- 
ence from the strong coupling theory is that at finite U 



there is a range of momentum q in which the lower edge 
of the two-holon continuum is described by two holons 
with the same velocity but different momenta. 

In order to investigate the behavior of the spectral 
weight of S{q,uj) near the edges of the spectrum, we 
relied on effective quantum impurity models. We have 
explicitly incorporated the 5*0(4) symmetry of the Hub- 
bard model at half filling by introducing SU{2) spinors 
for the high energy charge and spin modes. The internal 
degree of freedom in these spinors stems from degenerate 
particle and hole sub-bands. Once we have these objects, 
we write down effective Hamiltonians with marginal op- 
erators that are allowed by the spin and ry-spin SU{2) 
symmetries. In the effective impurity models the charge 
impurities are always decoupled from the low energy spin 
excitations due to symmetry. On the other hand, the spin 
impurities are coupled to the low energy spin excitations, 
but the coupling is marginally irrelevant due to Kondo- 
type physics. 

The operators that are associated with each threshold 
are also identified using symmetry. These operators must 
have the lowest scaling dimension that is allowed by the 
conditions that the excited state has the correct num- 
ber of impurities and that the operator has the correct 
quantum numbers for spin and 77-spin rotations. In the 
case of 5(q,cj), the operators are vectors of T^-spin and 
scalars of spin rotations. Due to the decoupling between 
low energy and high energy modes, the problem of edge 
singularities reduces to computing few-body propagators 
for the high energy part, which can be affected by final 
state interactions, and combining them with the correla- 
tion functions for the low energy part, which are known 
from conformal field theory. Simple power counting in 
the time decay of the total correlation function then de- 
termines the edge exponent /i for a given threshold. We 
have focused on S{q,u;), but the method can be read- 
ily applied to other dynamic response functions, such as 
the one-electron spectral function and the dynamic spin 
structure factor. 

The results of the effective quantum impurity models 
extend the validity of the low energy exponents'^^ = 1/2 
for q M and fi — 3/2 for 5 w tt to the regime of finite 
U , even though the spectrum is not relativistic as in the 
sine-Gordon model. The impurity models combined with 
the exact spectrum from the BA also provide the range 
of q over which these exponents hold. Remarkably, we 
found that the exponent /i = 1/2 at the lower edge of 
the two-holon continuum is verified only for q < qo{U), 
where qo(U) is determined by the infiection point of the 
holon dispersion relation. For q > q(^{U), there is a Van- 
Hove type square-root divergence along the lower edge 
of the two-holon continuum, due to the two holons that 
propagate with the same velocity but different momenta 
and do not scatter off each other in the integrable model. 
The existence of this divergent edge at finite U, near the 
zone boundary and at finite energies, is confirmed by the 
tDMRG results. Within the precision of the numerical 
results, we found no evidence for rounding of this sin- 
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gularity due to coupling to continuum below it, which 
would be apparent in the form of an exponential decay 
of the real-time correlation function. 

The agreement between the analytical predictions and 
the numerical line shapes obtained by tDMRG allowed 
us to explain how the line shape of S{q, uj) changes as a 
function of U, interpolating between the strong coupling 
and the weak coupling limits. Starting from strong cou- 
pling and decreasing U, we observed that the spectral 
weight inside the two-holon continuum decreases while 
the spectral weight below the lower edge of the two-holon 
continuum increases. The ?7 — > limit is nonperturba- 
tive, as expected from spin-charge separation and the 
Mott transition, and this is manifested in the dynamic 
response function through a discontinuous change in the 
edge exponents. For instance, while at J7 = S{q, ui) has 
a square-root divergence below the upper threshold of the 
electron-hole continuum, for arbitrarily small U this sin- 
gularity is removed and a square-root divergence forms 
above the lower threshold of the two-holon continuum. 

We end by commenting on the connection with exper- 
iments that show a sharp feature observed in the RIXS 
spectrum of ID Mott insulators for momentum near the 
zone boundary/^ This feature was interpreted as an exci- 
ton in Ref. 0, expected from the strong coupling theory 
for the extended Hubbard model, but broad two- 
holon resonance in Ref. 6. Our results for S{q,uj) of 
the integrable Hubbard model do not have any excitonic 
bound states, but also show a sharp feature near the zone 
boundary which is actually a square-root divergence at 
the lower edge of the two-holon continuum at finite U. 
Therefore, a possible interpretation of the experiments 
is that the sharp feature is the result of a slight round- 
ing of this Van Hove singularity in a system where the 
integrability breaking interactions (primarily the nearest 
neighbour interaction in the extended Hubbard model) 
are fairly weak. However, the nearest neighbor interac- 
tion is not guaranteed to be negligible since screening is 
typically rather weak in insulators such as Sr2 Cu03 . 
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Appendix A: Symmetry and elementary excitations 
in the Bethe ansatz solution 

Here we briefly discuss the relation of the operational 
representation of Ref. to the excitations considered in 
this paper. 

The pseudofermion dynamical theory^*' employs a uni- 
tary transformation originally devised to work in the 
strong coupling limit^^ that rotates electron operators to 
a basis where double occupancy is a good quantum num- 
ber. The rotated-electron configurations are then natu- 
rally expressed in terms of pseudoparticles whose discrete 
momentum values are BA exact quantum numbers. The 
occupancy configurations of the spin-1/2 spinous, rj-spin- 
1/2 r]-spinons, and spin-less and r^-spin-less c fermions of 
that representation generate both the representations of 
the spin SU(2) symmetry, 77-spin SU{2) symmetry, and 
charge hidden C/(l) symmetry algebras, respectively, and 
the model 4^ energy eigenstates. The spin-1/2 spinons 
are the spins carried by the rotated electrons of the singly 
occupied sites. The •q-spin-1/2 rj-spinons of projection 
— 1/2 and -1-1/2 refer to the Ty-spin degrees of freedom 
of the rotated-electron doubly occupied and unoccupied 
sites, respectively. The c fermions describe the charge 
hidden U{1) symmetry degrees of freedom of the rotated 
electrons of the singly occupied sites. The c fermion holes 
describe the hidden t/(l) symmetry degrees of freedom 
of the rotated-electron doubly occupied and unoccupied 
sites. 

The occupancy configurations of the spin-neutral com- 
posite sv fermions, each containing 2i> bound spinons, 
considered in Ref. 0, were called distributions of 
magnon bound states by M. Takahashi?^ Furthermore, 
the occupancy configurations of the 77-spin neutral com- 
posite rjv fermions of Ref. |4^ each containing 2iy 
anti-bound 77-spinons, correspond to his distributions of 
bound states of pairs. Specifically, the momentum oc- 
cupancy configurations of the c fermions, ?7-spin-neutral 
2z^-77-spinon composite 771^ fermions, and spin-neutral 2iy- 
spinon composite si fermions where v = l,...,oo gen- 
erate excitations described by the BA thermodynamic 
equations (2.12a), (2.12b), and (2.12c) of Ref. IH^, re- 
spectively. In units of 2'K/Na, the momentum values of 
those objects are the BA quantum numbers /j, J^"; ^''^^ 
in such equations, respectively. Here within the Ref4^ 
notation, the index n = v = 1, cx) in and refers 
to the number of anti-bound-?7-spinon pairs and bound- 
spinon pairs, respectively, and a = j is the momentum 
value index. 

Note that the two sets of BA thermodynamic equa- 
tions given in Eqs. (2.12b), and (2.12c) of Ref. 113, 
which are associated with 77-spin-singlet and spin-singlet 
excitations, respectively, have exactly the same struc- 
ture. This is consistent with the excitations described 
by the BA thermodynamic equation (2.12a) of that ref- 
erence referring to a degree of freedom other than rj- 
spin and spin. Consistently, in Ref. |4^ it is confirmed 
that the latter excitations generate representations of 
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the hidden C/(l) symmetry in the model extended global signs), we obtain 
[SO{A) X U{l)]/Z2 = [SU{2) X SU{2) x C/(l)]/Z| sym- 



metry. 

For the problem studied in this paper, only excita- 
tions generated by c momentum band and spin-neutral 
11 = 1 two-spinon si fermion band occupancy configu- 
rations play an active role. Those excitations also con- 
tain two 77-spinons, whose occupancies generate the three 
??-spin-triplet states. The spin-singlet excitations gener- 
ated by the two-spinon si fermion momentum occupancy 
configurations are described by the BA thermodynamic 
equations (2.12b) of Ref. [53 for n — 1 spinon pairs. 

In this paper we call holons and spinous the holes of 
the c fermion and si fermion momentum bands, respec- 
tively. Hence the spinous considered here are spin-neutral 
objects. This is in contrast to those of Ref. |4^, which 
carry spin 1/2. 

In the thermodynamic limit holons and spinous have 
dispersion relations edpc) and ea{ps)i respectively, where 
the dressed momenta pc,s and dressed energies ec,s are 
given by. 



-ps ■ 

(Al) 

The energy bands e^{q) and e^{p) and corresponding mo- 
menta q and p are given in Eqs. (A1)-(A4) of Ref. 0. 
For the present half filling case, the relation to Bessel 
functions provided in Eq. (A8) of that reference applies. 



Appendix B: Marginal coupling between chiral spin 
currents and spin impurity 

Consider the marginal operator in Eq. ((56)) . In this ap- 
pendix we derive the renormalization group (RG) equa- 
tions for this perturbation to the free Hamiltonian in Eq. 
(|5ip . The RG with high energy impurity modes is not 
standard, but the meaning is to investigate the effects of 
the perturbation when we approach the threshold where 
Hamiltonian (j5ip predicts a power-law singularity. The 
intuitive picture is that, as we approach the threshold, 
the energy of particle-hole excitations that the mobile 
impurity is allowed to scatter is reduced. Therefore, we 
shall consider the renormalization of the coupling con- 
stants Kfl.L when we integrate out an energy shell in the 
sub-bands near the Fermi surface. For consistency, the 
band width of the impurity modes must be reduced as 
well, but this effect will not be crucial for our conclusions. 

Let us focus on ki^ (the calculation for kr is completely 
analogous). We apply the perturbative RG.^^ The parti- 
tion function has the form 



Z = Tt exp 



-J d^x{n 



sn 



RLD) 



(Bl) 



Expanding for small (and omitting normal ordering 



l + 27rw,KL j d^x3L{x) ■ dI{x)tD^{x) 



{2ttVsKl) 



d'x / d'x' iT-),^,iT')u, 



(B2) 



X Jl {x)Jl {x')Dl, {x)D,^^^ {x)Dl, {x')D,^^ {x'] 



where Zq is the free part associated with Hamiltonian 
([51]) . The 0{hi\) term can generate corrections to kl 
when we integrate out "fast" modes. We use the operator 
product expansion of the spin currents^ 



I 

27rz 



e'^''V£(z') + 



(B3) 



where z = VsT -'rix is the complex argument of holomor- 
phic functions. In Eq. (jB2p . we must also take contrac- 
tion of Dg fields. For this purpose we need the impurity 
propagator in imaginary time 



(r.Z?,Ja:,r)i^I^.(0,0)) 



2^ 



*,j0(t)G(x,t), 



(B4) 



where K, is the momentum cutoff of the impurity sub- 
band. We obtain 

(B5) 

■nyUgT — ix) 

Note that we cannot take the limit /C 00 in Eq. (|B5p 
yet. (For the propagator in real time, this is possible and 
yields the delta function in Eq. ([5^ .) 

Using Eqs. (|B3| and HU in Eq. JHSJ, we find (keep- 
ing only corrections to k^) 



1 + 2ttVsHl I d^x3L -DlrD 



Tr{vsKL)^ / d^xJL-DlrD, 



Vgf + ix 



(B6) 



where {f,x) = (t — t' ,x — x') are the relative coordinates 
of the two points in Euclidean space-time. Importantly, 
the impurity propagates with a different velocity than the 
bosonic modes, thus the problem is not Lorentz invariant. 
Physically, this is more like a boundary problem, with a 
"mobile boundary" represented by the impurity that the 
bosonic modes have to track. Therefore, instead of a rota- 
tionally symmetric energy-momentum shell, we integrate 
out the "fast" modes contained in the strip — 00 < i < 00, 
1/A < |f I < 1/A', with A and A' = A - dA being the 
original and reduced energy cutoffs, respectively. The 
integration over x gives 



dx^^ 
Vsf + ix 



.1 



-K{Vs+Us)\f\ 



{Vs + Us)f 



(B7) 
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We can take the limit /C — >■ cxd in Eq. (jB7p . Moreover, 
we are integrating out short time differences f ^ 1 /A ^ 
I/Es, thus we can approximate e~^^'^ w 1. We are left 
with the imaginary time integral 



df 



2d£ 



(B8) 



l/A 



where d£ = dA/A. 

Finally, substituting the result in Eq. (jB6|) and reex- 
ponentiating, we find the RG equation for kl: 



(B9) 



dKL 
~df 



The RG equation for kji is obtained from Eq. (jB9P by 
the substitution kl — >■ hr,Us — > — Us. Since Ms < Ws, we 
conclude that k/, > and kji > are marginally irrele- 
vant. We believe this to be the correct sign for the cou- 
pling constants of the Hubbard model. Furthermore, we 
expect the marginally irrelevant hl,r operators to give 
rise to logarithmic corrections to edge singularities for 
SU{2) symmetric models, similarly to the effect in equal- 
time correlation functions. Logarithmic corrections are 
known to exist at the lower edge of the two-spinon con- 
tribution to the spin DSF for the Heisenberg model;^ 
but we do not pursue that calculation here. 



Appendix C: Exponent for threshold with two 
charge impurities and one spin impurity 

In this appendix we detail the calculation of the ex- 
ponent for the threshold U!2c2s{q<^ < Q < <Zo) in S{q,uj), 
which is described by two high energy holons and one 
high energy spinon, all moving with the same velocity. 
Other exponents can be obtained by similar methods. 

We find it convenient to calculate the exponent of 
S{q, uj) using the analytical continuation of imaginary 
time propagators to real time prescribed as follows. The 
zero temperature limit of the imaginary time propagator 
of the density operator is 



K(r)n_,(0))=^|(G5| 



--Egs)t 



(CI) 



Using the analytic continuation with the prescription 
IT — >■ — (t — ir]), we obtain 



|2„-i(S„-BGs)(*-i'7) 



(C2) 



where — >■ O"*" at the end guarantees the convergence of 
the sum. Taking the Fourier transform, we get 



+ 00 



dt e'^^^Ugit - ir])n^g{0)) 
27rY,\{GS\ng\iy)\H{uj^E, + EGs), (C3) 



which is the correct expression for S{q,u!). 

As argued in Sec. IIIIl due to yy-spin SU{2) symme- 
try the exponent for rj^ — excitation is the same as 
the exponent for the rj'^ — +1 excitation. In the lat- 
ter case we can treat the two holons as identical spinless 
fermions that do not interact via s-wave scattering. The 
only interaction in this three-body problem is between 
the holons and the spinon. In first quantization, we write 
down the effective Hamiltonian (for energies measured 
from the threshold) 



H^b — 



-vl 



2mc 

+ Vcs[S{xi 



~ X3) + S{X2 



P2 +P3) 
-2^3)], 



(C4) 



where particles 1 and 2 are the two holons and parti- 
cle 3 is the spinon, with canonically conjugated variables 
[xmPm] — i5nm- For a generic spinon-holon interaction 
potential, the parameter Vcs is related to the s-wave scat- 
tering length. The wave functions in the physical Hilbert 
space must be anti-symmetric with respect to exchanging 
1 and 2. The three-body propagator in imaginary time 
can be calculated from 



Gzb{x,T) = ($|e 



(C5) 



where P = pi + P2 + Pa is the total momentum operator 
and 



1$) 



Xl 



(8)|x3 = 0) 



:,X2 



Xl 



,X2 



(C6) 



is the initial state created by applying dl{x + §)rfj(a; — 
|)(i|(0) on the ground state. 

We perform a change of variables from xi,X2,x^ to 
X = [mc{xi+X2) + msX3,]/{2mc + ms),Xr = xi-X2,z = 
X1+X2 — 2x3 s-iid the associated conjugate momenta. The 
Hamiltonian becomes 



H: 



p2 



3b 



uP 



Pr 



Pz 



2(2mc + rus) 2mr 2mz 

+ 2Vcs[SiXr -Z)+ S{Xr + z)], 



(C7) 



where = mc/2 and ruz 



°.„^ I — In terms 

2{2ms+mc) 

of these new variables, the initial state has X = z = 
0, Xr — e. The x dependence of G3b{x) is entirely in 
the free centre-of-mass "particle". We note that rrir 7^ 
niz yinc,ms G IR. While > for holons below the 
inflection point, we assume > as well, which is 
easily verified in the strong coupling limit. 

First consider the simpler case Vcs — 0. In this case, 
all three particles are free and the propagator factorizes 



G3b{x,T) = Gcm{x,T)Gr{T)Gz{T). 



(C8) 



For the propagator of the centre of mass particle, which 
moves with velocity u, we shall use as in Eq. (|B5[) 



G cm {X-, t) 



sinh[/C(itT — ix)] 
Tr{uT — ix) 



(C9) 
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with cutoff JC <Si \{2'mc + rns)u\. For the other "particles", 
we have 



Gr{T) - / dkr sin2(fc^£)e-'='^/2'"- - 
for r ^ |m^|£^ and 

/OO 
-OO ^ T 



3/2 



(CIO) 



1/2 



Notice that Grir) decays faster because of fermionic 
statistics, which imposes that the wave function is an 
odd function of Xr- This is equivalent to the vanishing 
matrix element in Ref. 33. 

At the threshold W2c2s('?<i < Q < Q>), the three body 
propagator has to be combined with the low energy prop- 
agator of the chiral spinor, which has scaling dimension 
1/4. The integral over x gives (for Vs > u, there are two 
separate contributions from a pole and a branch cut in 
the lower half plane) 



dx 



Gcm{x,T) 1 

{vgT — ix)^/'^ [(us — u)r]^/^ ' 



in which we took the limit /C — )■ oo after the integration. 
Combining with Grlr) and G'z(t) and switching to real 
time IT ^ ~it ~ 'i-v) explained above, the remaining 
time integral gives 



where k is related to the energy hy E = fc^/2mr. 

Equation (jClSp describes the motion of a particle in 
two dimensions which is scattered by delta function po- 
tentials located along the lines Xr — ±2 cot a. We 
can solve the wave functions in the four regions of the 
{z,Xr) plane separated by these lines and then match 
the wave functions with a discontinuity in d^/d(j) at the 
boundaries. The solutions are of the form <E'(p, (j)) = 
(Cll) Er=o[^»e""^ -I- B„e-"'^]J„(fcp), where J„ denotes the 
Bessel function of the first kind. Imposing that the wave 
function is continuous everywhere and is anti-symmetric 
with respect to exchanging the two holons implies that 
<i>(/0 = 0, 0) = 0, hence = Bo = in all regions. 

In the long wavelength limit, /c — > 0, the delta func- 
tion potentials become impenetrable and the wave func- 
tion vanishes along the lines Xr = izcota. Importantly, 
these lines do not coincide with the z, Xf axis (in which 
the kinetic energy is diagonal) since rrir ^ m^. But we 
are interested in the behavior of the wave function for 
z = Q, Xr — e ^ i.e. approaching the origin along the 
The wave function already vanishes at p = 
due to the anti-symmetrization, therefore it is not af- 
fected by the delta function potential at a:^ = z = 0. 
As a result, for fee — >■ the eigenfunctions vanish as 
^{p = £,</> = 7r/2) ^ Ji(fce) ~ ke. This is the same 
behavior as obtained for Vcs = and leads to the expo- 
nent /X = 3/2 in S{q,uj) as in Eq. ((CT3)) . 



(C12) 



S{q,uj) 



iSujt 



dt 



{t - iTifl"^ 



^^)bu?l^. (C13) 



Now consider V^s 7^ 0. In this case the Xr and z parti- 
cles are scattered by the potentials in Eq. (jC7p . Nonethe- 
less, we argue that the edge exponent is the same as for 
Vcs = 0. First, we note that the exponent depends on 
the long-time behavior of G3b{x,t), which in turn de- 
pends on the behavior of low energy eigenfunctions for 
z = 0, Xr = e ^ 0. The extra power of 1/r in Eq. (jClOP 
is a result of the wave function vanishing as ^ fc^e for 
krS — >■ 0. Then we must ask whether Vcs 7^ modifies 
the behavior of the wave function in the long wavelength 
limit. 

Rescaling z — >■ ztana, — ^ p^cota with tana = 
y/rrir/mz 7^ 1 in Eq. (jC7p . the Xr and z part of the 
Hamiltonian becomes 

H2b = ^^^-ti^ + 2VcsSixr ± z tan a). (C14) 
2mT- 

We then introduce polar coordinates z = p cos 0, Xr = 
psmcj). The two-dimensional Schrodinger equation for 
the wave function $(0;^, z) = $(/0, </>) reads 



I d f 
p dp V dp 
2mrVcs cos a 



1 

p2 



[d{(l) ±a) + d{(l)±a + tt)], (C15) 



Appendix D: Absence of scattering between 
distinguishable charge impurities in the integrable 
model 



In this appendix we show that the coupling constants 
V^ and V^ in Eq. ((63|) are fine tuned to zero as a result 
of the integrability of the Hubbard model. Here integra- 
bility is understood as the existence of an infinite num- 
ber of local conserved quantities in the thermodynamic 
limit. The simplest nontrivial conserved quantity of the 
Hubbard model was discovered by Shastry^^ and can be 
written as^ 



Q3 = ^ [(ic]+i^^Cj-l,^ + h.c.) 



, (Dl) 



where J^j^„ = iCj^^ a'^j-<y + ^^^^ current density op- 

erator for electrons with spin a. The conserved quan- 
tity is almost equal to the energy current operator, 
differing only by a factor of 2 in front of C/<^ The en- 
ergy current operator — Y^- is defined from the 
continuity equation of the Hamiltonian density. Writing 
H ~ hj with 

h, = -(c]c^.+i+h.c.) + C/ (^n^-^ - ^ (^n,,; - ^ , (D2) 
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we obtain by taking the commutator of hj with H, 
which has the form of a discretized divergence 



The operator can be written as 

Q3 = 2j7^ + r, 



(D3) 



(D4) 



where Y = X), ^" = Z]j(~*4+iS-i + ^-C-)- Interest- 
igly, y is independent of \J and its density appears 
in the the commutator of the charge current density 
Jj — ^r, Jj,a with the total charge current J = ^jJj- 



ms 



(D5) 



We want to impose the conservation of in the ef- 
fective model Eq. A similar idea has been applied 
to the XXZ model)^ in which case it was shown that 
conservation laws lead to constraints on irrelevant op- 
erators at low energies, with consequences for dynamic 
correlation functions. Since the impurity model is phe- 
nomenological, we need a prescription to construct the 
conserved quantity directly in the field theory. The key 
is to use the continuity equations and relations (|D4|) and 
(jPSP since currents can be easily identified in the field 
theory. A caveat in applying Eq. (|D5P in the field the- 
ory is that the dimensions of the density of Y and 
differ by a factor of lattice spacing squared. This entails 
that when combining Y from Eq. inil with from 
Eq. (jD3P we must restore nonuniversal factors of short 
distance cutoff for dimensional analysis. 

The calculation of and Y in the field theory can 
be simplified using the local SU{2) algebra of DItD^ 
and D\tD^. The charge current density obtained from 
the continuity equation for the charge density n{x) 
Dl{x)T^-D,{x)+bl{x)T^-D,{x) is 

J{x) = dI{x)t' (uc~—d.^D^{x) 
\ rric J 

+DI{x)t' (u, - -^d,) D,{x). (D6) 

The commutator of the charge current density with the 
integrated charge current J — j dx' J{x') gives 



Dl{x) — - —d^ D,{x) 



nic mi: 



+ Dl{x)(^-^d^]D^{x) 
rrir mi 



(D7) 



Comparing with Eq. (IDSp . we conclude that the contin- 
uum version of y is y = J dx Y{x) with density 



Yix) = -Dlix) ( ^ - ) DAx) 

-Dlix) ( ^ - ) (D8) 

mr mi 



The energy current operator is obtained from the com- 
mutator iinlx),! dx"H{x')] = d^J^{x). We find 



cDI \Uc da 

mr 



Dr 



^a. D 



+ecD 

+2u,VfDlTD^-DlTD^, 



(D9) 



where we neglect operators with dimension higher than 
2. 

Using Eq. (jD4p , we construct the density of the con- 
served quantity 



i{x) 



2ErDl 



-2ErDl 



mr 



+AurV^DlDMD, 

+AUrVr^DlTDr ■ DItD, 
,.2 A 



(DIG) 



where Ec — ec~ l/{2mca''^), Ec — ic— l/(2mca^), with 
a the short distance cutoff. The density of Q3 in Eq. 
(jPlOp contains all the operators up to dimension 2 that 
are invariant under Ty-spin rotation but with different co- 
efficients than the Hamiltonian (l63l) . In fact, has the 
same symmetries as the Hamiltonian except for the sig- 
nature under parity transformation (parity symmetry is 
broken by hand in the effective impurity model by the 
definition of the impurity sub-bands). 

Finally, taking the commutator of Q3 = J dx Q3 (x) 
with H = J dx 'H{x), we are left with two dimension- 
three operators that do not not vanish in general 



\ mr mr 



C 



dxd^{DlDMDr 



+ ^^1 dxd,{DlTD,)-blTb, 



(Dll) 



We note that other terms cancel because the two sub- 
bands have the same velocity u,.. Recall that along the 
boundary line 0J2c{l > ^o) we have mc > and fhc < 
since the sub-band is below the inflection point of 
the holon dispersion and the Dc sub-band is above it. 
Moreover, > |toc| because the curvature of the holon 
dispersion is smaller close to the band minimum. Thus 
we have 



mr rhr 



mr 



1 



1 



2\mr\cP- 2mrCp' 



>0. (D12) 



The only way to ensure that the commutator in Eq. 
(|DTT|) vanishes is to set = Vf = 0. Therefore, the 
existence of a conserved quantity represented in the field 
theory by an operator of the form in Eq. (jPlOp requires 
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that there is no scattering between Dc and Dc holons. 
Importantly, integrabihty does not have any imphcations 
for the interaction between two impurities within the 
same spinor {Vg in Eq. and Vc in Eq. ([CTI) ). This 

fohows from taking ic = £c and rhc = rric in Eq. (jPlip . 
in which case the commutator vanishes identicahy. 

We also remark that the effective model in princi- 
ple also contains irrelevant interactions that have the 



same dimension (three) as the parabolic dispersion term. 
These irrelevant interactions, which were omitted in Eq. 
([63| . can contribute to the coefficient of the last two 
terms in the conserved quantity in Eq. (jPlOp . However, 
such terms do not contribute to the commutator in Eq. 
(jPlip (at the level of dimension-three operators), thus 
our conclusion is not affected by irrelevant interactions. 
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